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ABSTRACT 

A surprising finding of recent studies is the large number of Active Galactic Nuclei (AGN) associated with moderately massive 
black holes (log(M,/Mo) ~ 6 — 8), in the first billion years after the big bang (z > 5). In this context, a relevant finding has 
been the large fraction of candidate dual AGN, both at large separations (several kpc) and in close pairs (less than a kpc), likely 
in the process of merging. Frequent black hole merging may be a route for black hole growth in the early universe; however, 
previous findings are still tentative and indirect. We present JWST/NIRSpec-IFU observations of a galaxy at z = 7.15 in which 
we find evidence for a log(M./Mo) ~ 7.7 accreting black hole, as traced by a broad component of H 6 emission, associated 
with the Broad Line Region (BLR) around the black hole. This BLR is offset by 620 pc in projection from the centroid of strong 
rest-frame optical emission, with a velocity offset of ~40 km s^. The latter region is also characterized by (narrow) nebular 
emission features typical of AGN, hence also likely hosting another accreting black hole, although obscured (Type 2, narrow-line 
AGN). We exclude that the offset BLR is associated with Supernovae or massive stars, and we interpret these results as two 
black holes in the process of merging. This finding may be relevant for estimates of the rate and properties of gravitational-wave 
signals from the early universe that will be detected by future observatories like LISA. 


Key words: galaxies: high-redshift — galaxies: active — galaxies: supermassive black holes. 


1 INTRODUCTION 


After one year of operations, one of the most surprising findings in the 
field of extragalactic astronomy produced by the JWST (Gardner et al. 
2023; Rigby et al. 2023) is the unexpectedly high abundance of active 
* E-mail: hu215 9 cam.ac.uk galactic nuclei (AGN) during the first billion years of cosmic history 
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(z 2 5.5; Greene et al. 2023; Harikane et al. 2023; Kocevski et al. 
2023; Matthee et al. 2023; Yang et al. 2023; Maiolino et al. 2023). 
These are associated with accreting black holes and are identified 
through the properties of the spectral energy distribution (SED; e.g. 
Yang et al. 2023), through the presence of (narrow) nebular line 
emission typical of AGN (e.g. Scholtz et al. 2023; Übler et al. 
2023; Maiolino et al. 2024), and through the detection of a broad 
component of the Balmer lines, associated with the so-called Broad 
Line Region (BLR) of AGN (e.g. Greene et al. 2023; Harikane et al. 
2023; Kocevski et al. 2023; Matthee et al. 2023; Übler et al. 2023; 
Maiolino et al. 2023). With the strongest lines of the hydrogen 
Balmer series, Ho and H f, being redshifted out of atmospheric 
windows for z > 4 galaxies, many of these accreting black holes are 
now identified for the first time with JWST spectroscopy (primarily 
NIRSpec-MSA and NIRCam-slitless; Ferruit et al. 2022; Jakobsen 
et al. 2022; Boker et al. 2023; Rieke et al. 2023) — not least due to a 
remarkable lack of accompanying X-ray emission, previously taken 
as the main indicator to constrain the cosmic black hole accretion 
rate density (Madau & Dickinson 2014; but see also Bogdán et al. 
2023; Goulding et al. 2023). 

The JWST finding of the many AGN associated with moderately 
massive black holes has sparked the exploration of a large number of 
models, simulations, and theories in order to explain their properties 
and abundances (e.g. Mayer et al. 2023; Pacucci et al. 2023; Schnei- 
der et al. 2023; Trinca et al. 2023; Volonteri, Habouzit & Colpi 2023; 
Weller et al. 2023; Zhang et al. 2023; Di Matteo, Angles-Alcazar & 
Shankar 2023a; Bennett et al. 2024). These models involve various 
potential early black hole seeding scenarios, including: heavy seeds 
resulting from the direct collapse of primordial clouds; remnants 
of early massive stars; rapid merging of stars and (stellar mass) 
black holes in dense nuclear clusters. Episodes of super-Eddington 
accretion have also been suggested (Schneider et al. 2023; Trinca 
et al. 2023; Bennett et al. 2024). 

The role of merging between intermediate-mass black holes 
(IMBHs) has also been extensively investigated by models and 
simulations (Barai et al. 2018; Barausse et al. 2020; Volonteri et al. 
2020; Di Mascia et al. 2021; Valentini, Gallerani & Ferrara 2021; 
Mannerkoski et al. 2022; Vito et al. 2022; Volonteri et al. 2022; Chen 
et al. 2023; Di Matteo et al. 2023a, b). This is not only a potential 
channel for the growth of black holes in the remote universe, but also 
relevant for the detection of gravitational-wave signals from early 
cosmic epochs with future gravitational waves experiments. Within 
this context, significant progress has been made observationally in 
recent years through the detection of dual quasars and AGN separated 
by a few/several kpc at high redshift (z~1—3) (Mannucci et al. 2022; 
Ciurlo et al. 2023; Perna et al. 2023a), in a number significantly larger 
than predicted by cosmological simulations (e.g. De Rosa et al. 2019; 
Rosas-Guevara et al. 2019; Volonteri et al. 2022; Chen et al. 2023). 
At even higher redshift (z>4), JWST multi-object spectroscopic data 
have suggested the presence of a significant population of dual AGN, 
on scales smaller than 1 kpc, based on the complex profile of their 
BLR, suggesting the presence of two separate BLRs associated 
with two accreting black holes (Maiolino et al. 2023). However, 
in the latter case no spatial information was available to confirm the 
presence of two AGN. Moreover, the complex profile of the broad 
lines may potentially also be associated with some peculiar properties 
of an individual BLR. 

In this paper, we report the detection of a broad H£ line in 
the z — 7.15 galaxy system ZS7, which identifies it as a (Type 
1) AGN. ZS7 is also known by its CANDELS ID (Grogin et al. 
2011; Koekemoer et al. 2011; Nayyeri et al. 2017) COSMOS13679 
(Pentericci et al. 2016, 2018), COSY-0237620370 (Roberts-Borsani 


MNRAS 531, 355-365 (2024) 


et al. 2016; Laporte et al. 2017; Witten, Laporte & Katz 2023b), and 
COS-zs7-1 (Stark et al. 2017). This galaxy was already identified to 
host an AGN, based on the detection of N vA1239 and He1.1640 
in the rest-frame UV (Laporte et al. 2017). However, the peculiarity 
of this system is that the BLR is offset from the primary (narrow) 
nebular line emitter. Additionally, also the latter emitter displays 
nebular line signatures typical of an AGN. We interpret our results 
in the context of supermassive black hole mergers, which will be 
detectable with future gravitational waves observatories like LISA. 

Throughout this work, we use the AB magnitude system and 
assume a flat ACDM cosmology with Qm = 0.315 and Hy = 
67.4kms'! Mpc! (Planck Collaboration et al. 2020). With this 
cosmology, 1” corresponds to a transverse distance of 5.27 proper 
kpc at z = 7.15. 


2 NIRSPEC-IFU OBSERVATIONS AND DATA 
PROCESSING 


Raw data files were downloaded from the Barbara A. Mikulski 
Archive for Space Telescopes (MAST) and subsequently pro- 
cessed with the JWST Science Calibration pipeline! version 1.8.2 
under the Calibration Reference Data System (CRDS) context 
jwst-1068.pmap. We made several modifications to the default 
reduction steps to increase data quality, which are described in detail 
by Perna et al. (2023a) and which we briefly summarize here. Count- 
rate frames were corrected for 1/f noise through a polynomial fit. 
During calibration in Stage 2, we removed regions affected by failed 
open MSA shutters. We also removed regions with strong cosmic 
ray residuals in several exposures. Remaining outliers were flagged 
in individual exposures using an algorithm similar to LACOSMIC (van 
Dokkum 2001): we calculated the derivative of the count-rate maps 
along the dispersion direction, normalized it by the local flux (or 
by three times the rms noise, whichever was highest), and rejected 
the 96" percentile of the resulting distribution (see D'Eugenio et al. 
2023, for details). The final cube was combined using the ‘drizzle’ 
method, for which we used an official patch to correct a known bug.? 
The main analysis in this paper is based on the combined cube with 
a pixel scale of 0.05". We used spaxels away from the central source 
and free of emission features to perform a background subtraction. 

ZS7 was observed with the NIRSpec instrument (Jakobsen et al. 
2022) in Integral Field Unit (IFU) mode (Bóker et al. 2022) onboard 
JWST as part of the NIRSpec-IFU GTO programme ‘Galaxy Assem- 
bly with NIRSpec IFS’ (GA-NIFS) under programme 1217 (PI: Nora 
Lützgendorf). The NIRSpec data were taken on 2023 May 9, with a 
medium cycling pattern of eight dithers and a total integration time of 
4.1 h with the medium-resolution grating/filter pair G395M/F290LP, 
covering the wavelength range 2.87 — 5.27 um (spectral resolution 
R ~ 700 — 1300; Jakobsen et al. 2022). It was also imaged with the 
NIRCam instrument (Rieke et al. 2023) in eight filters though the 
PRIMER programme (PID: 1837, PI: James Dunlop). 


3 DATA ANALYSIS 


3.1 Spectral fitting and detection of broad H £ 


Fig. 1 shows two spectra (zoomed around H f) extracted from 
the IFU cube at the location of the [O m]à5007 centroid and 
at a location about 0.12" to the north-west, where a broad 


'https://jwst-pipeline.readthedocs.io/en/stable/jwst/introduction.html 
*https://github.com/spacetelescope/jwst/pull/7306 
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Figure 1. Top left: Colour-composite image of ZS7 and its immediate environment created using public NIRCam data from the PRIMER programme (PI: 
James Dunlop). The approximate NIRCam PSF is indicated by the white filled (F200W, ~0.06”) and open (F444W, ~0.14”) circles. The effective NIRSpec-IFU 
field of view (FOV) is indicated by the dashed white box. Top right: Zoom-in on the environment of ZS7, as indicated by the rectangular solid white box in the 
top left panel. The colour scale shows a line map integrated over 90 A around the [O 11145007 position, where we indicate the positions of two faint companions 
detected spectrally at the same redshift with dashed grey circles. The centroids of the [O 11]A5007 emission and the BLR emission are indicated by a white 
plus and a golden star, respectively, where the small white circles on top of the symbols indicate five times the fitting errors on the centroid positions. The 
approximate NIRSpec PSF at 4 jum is indicated by the white circle. The centroid of the [O 11]A5007 emission coincides with the NIRCam F200W emission, 
while the BLR location coincides with the F444W emission, reflected in the colour gradient in the NIRCam false-colour image. Bottom: Zoom-in on the spectra 
extracted from a single spaxel at the [O 11]4.5007 centroid position (left) and at the BLR location (right). We show our full fit in blue, with narrow emission 
lines in green and the BLR component in purple. At the BLR location, the broad H B emission is evident, while the [O 11]AA4959, 5007 emission is narrow at 
all locations. Some flux of the BLR point-source is spread through the point-spread function (PSFrwym ~ 0.13 — 0.14" at 4 um) and therefore is still visible at 


the [O 111]45007 centroid position. 


component of H is clearly visible. We fit the spectra with a 
combination of Gaussian functions for the line emission, and a 
power-law continuum. We include a narrow component for each 
of the following lines: [O1JA3727, [O1]A3729, [Ne m]à3869, 
He 1A3889+H8, [Ne1J]4A3967, He, Hô, Hy, [O m]à4363, H f, 
[O 11]A4959, [O m]à5007. The intensity ratio of the [O m] doublet is 
fixed to [O 11]45007/[O 111]A4959 = 2.98, based on atomic physics 
(Storey & Zeippen 2000). We note that the [O rr] doublet as well as 
[Ne 11]43967 and He are blended. We show our full fits to spectra 
extracted from the BLR location and the centroid of the [O 111]4.5007 
emission in Appendix A. 

As mentioned, for H £, a broad component is observed in some 
spaxels. We see no evidence for broad components in other lines. As 
it will be discussed later on, we interpret this broad H 6 component 
as a BLR, offset relative to the [O 111]A5007 centroid. 

Following Übler et al. (2023), we rescale the formal noise based 
on the ‘ERR’ extension with a measurement of the standard deviation 
in regions free of line emission to take into account correlations due 


to the non-negligible size of the PSF relative to the spaxel size. This 
increases the noise by a factor of up to two for individual spaxels. 

We first fit the integrated spectrum extracted from 3 x 3 spaxels 
centred on the BLR location to robustly constrain the shape of the 
BLR emission. Measuring the BLR FWHM on a larger aperture 
including all spaxel with S/Ng g, grr = 5 (see Fig. 2) gives consistent 
results within the uncertainties. Here and elsewhere, the S/N is 
measured directly from the best fit. Then, we fit the emission 
extracted from individual spaxels by fixing the BLR FWHM to the 
previous best fit value. We list the resulting emission line fluxes 
extracted at the BLR location and [O r1]45007 centroid position in 
Table 1. Zoom-ins on the corresponding spectral fits extracted from 
individual spaxels are shown in Fig. 1. 

If we integrate over a larger aperture centred on the [O 11]4.5007 
emission in ZS7, we can detect additional fainter lines. A circular 
aperture with radius 0.3" (encompassing both the [O 111]45007 emitter 
and the BLR, similar to the region shown in Fig. 2) reveals tentative 
line emission at the expected wavelengths of H11, H10, H9, and 
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Figure 2. Top left: Continuum flux density map at A = 2.9 — 3 um. Bottom left: Continuum flux density map at A = 4.9 — 5 wm. Top centre, bottom centre, and 
bottom right: Integrated line flux maps with S/N > 5 from our spectral fitting for the narrow [O 11]45007 component, the H 8 BLR component, and the narrow 
[O 11]44363 component, respectively. Top right: Velocity map of the narrow [O 111]AA4959,5007 emission. In each panel, the black cross and star indicate the 
centroid of the [O m]à5007 emission and the BLR location as shown in Fig. 1, respectively. The black scale bar to the left indicates 0.2" ~ 1 kpc, and the grey 
outline indicates the mask used for our full spectral fitting. There is a clear shift in the location of the continuum emission from bluer to redder wavelengths, 
aligning with the centroids of the [O 111]4.5007 emission and the BLR position, reflecting the colour gradient seen in the NIRCam image (Fig. 1). The flux of the 
auroral [O 111]44363 line peaks at the BLR position, but covers also the [O 111]A5007 flux peak, indicating a clear gradient of the [O 11]4.5007/[O 111]44363 line 


ratio from the position of the AGN to the [O 111]A5007 emitter. 


Table 1. Observed narrow emission line fluxes extracted from one spaxel 
in units of 107% erg s~! cm’? from full spectral fits to the BLR location and 
[O 11]A5007 centroid position and derived quantities. 


BLR location [O 111]45007 centroid 
[O 1]AA3727, 3729 4.8 + 0.8 5.8 + 0.6 
[Ne m]à3869 4.6 + 0.5 7.1 + 03 
He I+H8 0.9 + 0.4 2.7 + 0.3 
[Ne 11]A3967+H € 2.6 + 0.7 3.9 + 0.5 
Hó 1.7 + 04 41+ 0.3 
Hy 2.9 + 0.4 5.5 + 0.3 
[O 11]A4363 3.1 + 0.4 2.6 + 0.3 
Hf 6.2 + 0.5 14.1 + 0.5 
[O m]à5007 72.9 + 1.8 117.1 + 1.5 
12+log(O/H) 7.69 + 0.08 7.81 + 0.04 
Te (0m]s5007 [K] 22700 + 1900 17400 + 500 


a detection of He 145876. While He 1141640 has been detected in 
previous observations from the ground (Laporte et al. 2017), we 
do not detect He 144686 in ZS7. This is not necessarily surprising, 
because He 1141640 is the brighter line. If we extract a spectrum 
collapsed from a pseudo-slit corresponding to the previous X-Shooter 
observations by Laporte et al. (2017), we find an upper limit on the 
He 144686 flux which is consistent within the uncertainties with their 
measurement for He 14.1640 within a factor of 8-13, consistent with 
theoretical expectations. 


3.2 Maps and velocity field 


We use the fluxes resulting from full spectral fitting to obtain emission 
line maps, for the regions where the lines are detected at more than 
5c. In Fig. 2 we show specifically the flux maps of [O 11]45007 
(which is also shown in Fig. 1), the broad Hf component, and 
[011]44363. 
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The map of the broad component of H f is consistent with the 
PSF, that is, a spatially unresolved source. The maps confirm what 
was already visible from Fig. 1, namely the H 8 broad component is 
clearly offset from the [O m]à5007 line centroid. 

The auroral line [O 11]44363 peaks at the location of the H 6 
broad, not unexpectedly, as strong [O 111]A4363 emission was already 
pointed out to be likely associated with AGN activity (Perna et al. 
2017; Dors et al. 2020; Brinchmann 2023). However, it is interesting 
to see that it also extend to the [O m]à5007 peak; we will discuss 
that the strength of the [O 111]A4363 line also at this location suggest 
the presence of an additional (Type 2) AGN. 

We also derive the [O 111]45007 velocity map, as shown in Fig. 2. 
We use the same fitting approach as described in Section 3.1, but 
fit only the [O11]A44959,5007 emission which has the highest 
signal-to-noise ratio. The velocity is derived from the line centroid 
after subtracting the systemic velocity measured at the [O 11]4.5007 
centroid position. From our fits, we find an average uncertainty of 
6km s^! on the velocities shown in Fig. 2. The region around the BLR 
location appears clearly separated in velocity space, as indicated by 
the cluster of positive relative velocities. The absolute difference 
in velocity to the centroid of the [O 11]A5007 emitter is only about 
20 km s^!, but significant. Being at the the same redshift, the small 
projected distance and the small velocity difference suggest that we 
are witnessing the imminent merger of two systems, one of which 
hosts a BLR AGN. 

Finally, we also extract and map the continuum in the bluest (2.9— 
3 um) and reddest (4.9-5 um) parts of the spectrum. As illustrated 
in Fig. 2, the red continuum peaks at the location of the broad 
H f, while the blue continuum peaks on the [O 11]A5007 centroid. 
This essentially confirms what was already seen in the NIRCam 
colours shown in Fig. 1. This further indicates that we are likely 
observing two distinct objects in close proximity. The blue continuum 
centered on the [O 111]A5007 centroid is likely tracing a young stellar 
population. Instead, it is not straightforward to interpret the nature of 
the compact red source at the location of the broad H £: it could be 
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a reddened or aged stellar population, but it may also be associated 
with the reddened continuum of the AGN accretion disc. It is not 
possible to easily distinguish between these two scenarios, because 
the continuum spectral information is limited to 3600 À res < 
5200 À and has modest S/N. 


3.3 Stellar mass estimate for the [O n1]à 5007 emitter 


We estimate the stellar mass of the [O 11145007 emitter through a full 
SED modelling using imaging data from the PRIMER programme 
on the pixel level, following Pérez-González et al. (2023) (see also 
D'Eugenio et al. 2023; García-Argumánez et al. 2023). SEDs were 
constructed for each one of the 26230 mas (150 pc, with the PSF 
FWHM of the reddest image being around 5 pixels) pixels of the 
system counting with more then three NIRCam filters with S/N > 
5. The NIRCam data covers the wavelength range 0.9-4.4 wm with 
eight filters (one medium-band filter, the rest being broad-band). 
The SEDs were fitted to stellar population models by Bruzual & 
Charlot (2003), assuming a delayed exponential for the star formation 
history (SFH) with time-scale t taking values from 1 Myr (i.e. an 
instantaneous burst) to 10 Gyr (i.e. a constant SFH). The attenuation 
was treated with the Calzetti et al. (2000) law, with A(V) values from 
0 to 5 mag. The metallicity was left to vary between 0.005Z_ and 
Zo. The IMF was set to Chabrier (2003). With a total of four free 
parameters, the best-fitting model (using a Monte Carlo method to 
estimate errors) is scaled to the photometry to get an estimation of 
the (surviving) stellar mass. 

To estimate the stellar mass of the [O 11]4.5007 emitter, we exclude 
the Northern region where the BLR is detected with S/N > 5 to 
mitigate contamination from the AGN (see Fig. 2). This leads to a 
total stellar mass of M, ~ 2.9 + 0.9 x 10? Mo. If we fill in the 
excluded region with values from the Southern part of the galaxy 
at similar distance from the [O 11]A5007 centroid, we find M, ~ 
3.8 + 1.1 x 10? Mo. This estimate is consistent with what has been 
found in previous work from SED fitting to the CANDELS and 
also IRAC data (Pentericci et al. 2016; Laporte et al. 2017). This 
shows that the [O 111]4.5007 emitter is a moderately massive galaxy 
for its redshift (see e.g. Curti et al. 2023; Nakajima et al. 2023), 
expected to host an intermediate-mass to massive black hole in its 
centre. 

In Section 7, we will argue that the narrow line emission associated 
with [O 11]44363 may also trace recent or ongoing AGN activity. 
Therefore, we make another estimate of the stellar mass of the 
[O 11]A5007 emitter by excluding not only the region where the 
Hf BLR is detected with S/N > 5, but also the region where the 
narrow [O m]à4363 is detected with S/N > 5. This would lead to 
a stellar mass of M, ~ 1.7 + 0.5 x 10? Mo (without filling in the 
excluded regions with other values). We consider this a lower limit 
on the total stellar mass of the [O 11]A5007 emitter, and note that 
a reduction of the stellar mass by a factor 2-3 does not impact the 
results and conclusion presented in the following sections. 


4 EVIDENCE FOR AN OFFSET TYPE 1 AGN 
AND ENVIRONMENT OF ZS7 


Ata projected distance of 620 pc from the centroid of the [O 111]A5007 
emission in ZS7, we detect a broad component at the wavelength 
of Hf (Fig. 1). As discussed, the emission is consistent with 
stemming from a point source, spread by the PSF (PSFrwum ~ 
0.13 — 0.14” at 4 jum; see Fig. 2). The broad component does not 
have a counterpart in the forbidden lines, in particular not in the 
strong [O1]AA4959, 5007 doublet, hence it cannot be associated 


with emission of high-velocity ISM or CGM, for instance due to 
galactic outflows. Therefore, the most plausible explanation is that 
the broad H £ is tracing the BLR of an accreting black hole, that 
is, a Type 1 AGN (we discuss and discard other interpretations in 
Section 6). 

The finding of an AGN in ZS7 is not so surprising, as evidence 
for an AGN was already found in the past based on the detection of 
N VA1239 and He 1141640 in the rest-frame UV (Laporte et al. 2017). 
The novelty of our finding is that the Type 1 AGN is offset relative to 
the [O n]à 5007 emitter. Also based on the kinematic offset around 
the BLR (Fig. 2 and Section 3.2) and the markedly different colours, 
these results suggests that the Type 1 AGN traced by the broad H f 
is hosted in a system separate from the galaxy traced by the bulk of 
the [O 111]A5007 emission, and in the process of merging with it. The 
ongoing merger may be channelling gas to the 5 x 107 Mo black 
hole as a consequence of tidal interactions (Di Matteo, Springel & 
Hernquist 2005; Johansson, Naab & Burkert 2009; see also e.g. 
Dougherty et al. 2024). 

We consider the possibility that the BLR is actually at the centre of 
a single system but partially obscured by dust, while the [O 11]A5007 
brighter emission region is tracing the less obscured part of the galaxy 
(or a less obscured part of an extended narrow-line region). However, 
the dust extinction towards nebular lines as traced by the H 8/H y 
Balmer decrement does not show a significant gradient between 
the two locations and, if anything, slightly lower reddening towards 
the narrow nebular lines at the BLR location (see Appendix B and 
Fig. B1). In addition, the fact the [O 111]A4363 line, which is clearly 
resolved across the two locations, peaks at the BLR location, provides 
further evidence that the BLR location is not the dust-obscured centre 
of the [O 111]45007 emitter. In the case of differential reddening, one 
would expect that [O 11]A4363 would peak at the same location as 
[O m]à 5007, whereas it peaks at the BLR location. 

It is also implausible that the BLR is tracing the centre of a 
single system and that the [O 11]A5007 centroid is tracing a burst 
of star formation in its disc, at ~600 pc from the nucleus. Indeed, the 
[O 11]45007 emission is not an unresolved clump of star formation, 
it is a well-resolved structure extending over ~1.5 kpc. Additionally, 
if the BLR was the centre of a single system then the velocity field 
should be approximately symmetric relative to it; on the contrary, 
the velocity field around the BLR stands out for having a clear offset 
relative to the surrounding region, just as expected for a separate 
system. Furthermore, there is no evidence of significant extended line 
or continuum emission north-west of the BLR. Finally, we will show 
that also the [O m]à5007 centroid is likely hosting another (Type 2) 
AGN, indicating that it is associated with a separate nucleus. 

At the same redshift as ZS7, we detect three fainter companions 
through emission lines of H £ and [O 11]A4.4959, 5007: two to the 
North at projected separations of 2.5 and 4.1 kpc, and one to the 
south-east at a projected distance of 3.9 kpc. Only the closest satellite 
is visible in the imaging data (see Fig. 1), and has been previously 
photometrically identified by Witten et al. (20232). The high number 
of companions, in addition to the offset BLR, suggests that we are 
witnessing a snapshot of the complex assembly of a larger system 
just 740 million years after the big bang. 


5 BLACK HOLE MASS AND LUMINOSITY OF 
THE OFFSET TYPE 1 AGN 


By assuming standard virial relations, we can estimate the black 
hole mass using our best fit values of the H 6 BLR flux and FWHM, 
based on the 20-spaxel integrated spectrum around the BLR location, 
where S/Nug BLR = 5. The FWHM is corrected for the instrumental 
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Table 2. Measurements of the broad Hf component extracted from an 
integrated spectrum with S / Nu ,BLR)5 (see bottom centre panel of Fig. 2), 
and corresponding black hole properties. 


Measurement Value 
FWHM broad [kms] 3660 + 300 
Fup, broad [erg s! cm ?] (1.7£0.2) x 10715 
Lpp,broad [erg s 1] (1.0 £ 0.1) x 10 
log(Mau/ Mo)Hg broad 71.7404 
log(Looi/Cerg/s))ng broad 449+0.1 
log(Lraa/(erg/s)) 45.8707 

ÀEdd 0.14 


spectral resolution which is approximately R ~ 1000 at the observed 
wavelength of H £, corresponding to o ~ 127 km s! (Jakobsen et al. 
2022). Using the single-epoch calibration by Greene & Ho (2005) 
derived from data of nearby AGN, we find a black hole mass of 
log (M,/Mo) = 7.7 + 0.4. The uncertainties are dominated by the 
scatter of local scaling relations, which we take into account by 
adding 0.4 dex in quadrature to our error budget (e.g. Ho & Kim 
2014). We also estimate the AGN bolometric luminosity to be Lo) = 
8.5(+ 0.8) x 10% erg s! based on the H £ BLR flux following Dalla 
Bontà et al. (2020), indicating that the black hole is accreting at about 
14 percent of its Eddington limit. We report these measurements in 
Table 2. 


6 EXCLUDING ALTERNATIVE SCENARIOS 
FOR THE BROAD H f: SUPERNOVAE, 
OUTFLOW, AND MASSIVE STARS 


While we recall that outflows as fast as a few 1000 km s! are only 
associated with AGN feedback (e.g. Fabian 2012; Heckman & Best 
2014; Carniani et al. 2015; Forster Schreiber et al. 2019; Veilleux 
et al. 2020; see also discussion by Maiolino et al. 2024), the outflow 
scenario should result in even more prominent broad emission in 
[O 11]AA4959, 5007, as observed in galactic outflows locally and at 
high redshift, and especially because outflows are typically metal 
enriched with high [O 11]A5007/H 6 (Origlia et al. 2004).? 

We briefly discuss alternative interpretations of the broad H £ 
component: a local outflow, a supernova (SN) in the galaxy disc, or 
a cluster of Wolf-Rayet or very massive stars. 

Core collapse SNe do show broad Balmer line emission; at 
their peak luminosity their broad Hf luminosity is of the order 
of ~10*° erg s^! , although some extremely luminous SNe can reach 
broad H f luminosities^ of a few times ~10*! ergs! (Pastorello 
et al. 2002; Fransson et al. 2014; Tartaglia et al. 2016; Kokubo et al. 
2019; Dickinson et al. 2023), still lower than what is observed in 
ZS7 (Lug, broad ^v 10? erg s; Table 2). The additional important 
difference is that when SNe reach such high luminosities, the broad 
component of H « and H £ is generally much larger than observed in 
ZS7, typically several 1000 km s^! and even exceeding 10 000 km s^. 
Moreover, core-collapse SNe, especially the luminous ones, have 
an asymmetric blueshifted profile (from a few 100 km s! to a few 
1000 km s!) and/or P-Cygni profiles (e.g. Gal-Yam 2012; Taddia 
et al. 2013; Gutiérrez et al. 2017), while the broad H 8 component 


3We note that in high-metallicity systems (12+log(O/H)=, 8.5) [O m]à5007 
can be suppressed with respect to H £ (e.g. Curti et al. 2020). However, our 
metallicity estimates for ZS7 are much lower (see Table 1). 

4We note that SN luminosities in the literature are generally reported for Ha, 
and H &/H £ can be quite large. 
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in ZS7 is symmetric (Fig. 1). Since we detect the continuum, we can 
also rule out strong absorption lines. We also note that the spectrum of 
core-collapse SNe (especially in the early phases) is further generally 
characterized by prominent [O 1]46300 emission and/or Fe bumps 
near H £, which are not seen in ZS7. 

Additional compelling evidence against the SN scenario comes 
from the fact that in this case the source at the location of the broad 
H £ should display some flux variation. We can leverage the fact 
that the NIRCam and NIRSpec-IFU observations were obtained at 
two different epochs separated by 5 months, that is, 20 d rest frame. 
While this is not a long-time lag, depending on the phase of the 
putative SN light curve, we would expect a flux variation ranging 
from a factor of several to at least 20 per cent. We have consistently 
extracted the flux from the same aperture both from the NIRCam 
and NIRSpec data. For NIRCam we have used the filters overlapping 
with our NIRSpec spectrum, that is, FA44W, F410M, and F356W, 
and in NIRSpec we have extracted the spectra matching the same 
wavelength range as the three filters. At least a 5 percent cross- 
calibration uncertainty is estimated due to the fact that the PSF of 
the two instruments are not identical (i.e. slight differential losses, 
despite the large adopted aperture), as well as slight background 
disuniformities in the IFU data. Nonetheless we obtain that the fluxes 
in the three bands are fully consistent within 3 per cent, 2 per cent, 
and 0.07 per cent, respectively. Such a good match, and the implied 
constancy of the flux (both of the continuum and of the broad H £, 
which is included in all three filters), makes the single-SN scenario 
simply impossible. 

However, one can explore whether the broad H f is not due to 
a single SN, but the overlap of multiple, colocated SNe associated 
with star formation within the small region associated with the BLR. 
However, even assuming the extreme case of no AGN contribution to 
the narrow H £, this would give an upper limit on the star formation 
rate of «6 Mg yr", that is, « 107! SN per year; even in the extreme 
case that the peak of the broad-line emission of each SN lasts for 
one year, this would imply an average broad Hf luminosity of 
10?? =% erg s™!, much lower than observed in ZS7. Also, one would 
expect the peak of the core-collapse SNe to pile-up at the location of 
the peak of star formation, that is, at the [O 11]A5007 centroid. 

Although it is not a common feature, broad Balmer lines have been 
associated with clusters of Wolf-Rayet stars in some individual cases 
(e.g. in NGC 4214, see Sargent & Filippenko 1991). More distinctive 
features of Wolf-Rayet stars are for instance the ‘blue bump’ around 
4600—4720 A, associated with broad He 144686, N 1144640, and/or 
C rI/IV44650, and the ‘red bump’ around 5800 A, associated with 
CIVA5808 (e.g. Sargent & Filippenko 1991; Terlevich et al. 1996; 
Schaerer et al. 1997; Schaerer, Contini & Kunth 1999; Brinchmann, 
Pettini & Charlot 2008; Monreal-Ibero et al. 2017). These spectral 
signatures are not seen in ZS7. We calculate an upper limit on 
the flux in the ‘blue bump’ region, and on He1144686, from our 
aperture where the H 6 BLR is detected with S/N > 5. From this, 
we find ‘blue bump’/H £ < 0.033 and He 144686/H B < 0.006. For 
comparison, Sargent & Filippenko (1991) find ‘blue bump’/H B ~ 
0.91 in the region where they also detect broad Ha, and values 
of He 1114686/H 8 ~ 0.01 — 0.1 are detected in nearby Wolf-Rayet 
galaxies (e.g. Vacca & Conti 1992; Thuan, Izotov & Lipovetsky 1996; 
Fernandes et al. 2004). To further constrain from the H £ flux and 
the upper limit on He 1144686 the possible contribution from Wolf— 
Rayet stars to the spectral features in ZS7, we follow Equation 15 by 
Vacca & Conti (1992) to calculate the ratio of the number of WNL 
stars, the dominant subtype of Wolf-Rayet stars, and the ‘equivalent’ 
number of O stars, and find 0.019. Considering the detection of 
He 114.1640 in ZS7, the two independent observations discussed by 
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Laporte et al. (2017) reveal line widths that are consistent with 
the spectral resolution of the respective instruments, while again 
generally broad emission would be expected for Wolf-Rayet stars 
(e.g. Chisholm et al. 2019; Leitherer, Lee & Faisst 2019; Vanzella 
et al. 2020). In summary, while we cannot exclude the presence of 
Wolf-Rayet stars, we can conclude that their potential contribution 
to any spectral features in ZS7 is minor. This specifically holds for 
the broad H f emission. 

Very massive stars are also known to show broad Balmer lines, but 
in these cases the Equivalent Width of the broad component of H 8 
is only a few À (Martins et al. 2020), while in ZS7 the Equivalent 
Width of the broad component of H £ is about 60 A. Furthermore, in 
the case of very massive stars the broad H £ is also accompanied by 
prominent and broad He 1114686 (Martins & Palacios 2022), not seen 
in ZS7. Finally, we note that an outflow driven by a large numbers 
of massive stars has been invoked as a possible origin of broad off- 
nuclear Ha emission without a counterpart in [O 111].4.4959,5007 in 
a nearby galaxy (Rodríguez Del Pino et al. 2019). However, in the 
latter case the width of the broad Ha is only ~800 — 1000 km s^!. 
Moreover, in the latter scenario, as well as all previous scenarios in 
which the broad H f is associated with clusters of young stars, one 
would expect the broad H £ to peak at the location of the peak of 
star formation, that is, at the location of the peak of (narrow) nebular 
emission, and not offset as in the case of ZS7. 

Given the several arguments outlined above, and also the fact 
that previous studies already found evidence for an AGN in ZS7 
(Laporte et al. 2017), we conclude that the most plausible explanation 
is that the broad H 8 component is tracing the BLR of a massive, 
accreting black hole. Furthermore, the H 8 BLR and the [O r11]4.5007 
distribution are likely associated with two different systems in the 
process of merging. This was already discussed in Section 4, and 
will be further discussed in the following two sections. 


7 INDICATION OF A SECOND, TYPE 2 AGN IN 
THE [O 111]35007 EMITTER 


As discussed by various authors, the detection of Type 2 AGN, that 
is, for which the BLR is obscured along our line of sight, and only the 
Narrow Line Region (NLR) is visible, is more challenging at high 
redshift. Indeed, the narrow line ratios of AGN at z > 4 are displaced 
from the local AGN locus in classical optical diagnostic diagrams, 
such as the BPT (Baldwin, Phillips & Terlevich 1981), involving 
[O 11]45007/H £ versus [N 11]6584/H a, and the diagrams proposed 
by Veilleux & Osterbrock (1987), in which [N 116584 is replaced 
by the [S11]6720 doublet (Harikane et al. 2023; Kocevski et al. 
2023; Scholtz et al. 2023; Übler et al. 2023; Maiolino et al. 2023). 
This is likely a consequence of the gas in the NLR becoming metal 
poor, which results in an overlap of AGN with the sequence of star- 
forming galaxies on these same diagrams. Some of the alternative 
diagnostic diagrams to identify AGN based on UV lines proposed 
by, for instance, Scholtz et al. (2023), are not usable in our case, due 
the limited spectral coverage of our data. 

However, an additional diagnostic can be provided by the strong 
[O 11144363 auroral line, which is found to peak at the same location 
of the broad H £ (Fig. 2), with an intensity comparable to H y. Such 
high [O 111]A4363 is very rare and, in combination with [O m]à 5007, 
indicates very high-gas temperatures, of the order of T; fo mjas007 = 
22700 + 1900 K (following Dors et al. 2020). Such a high ratio of 
[O 11]44363/H y is not seen in nearby star-forming galaxies, but it 
is seen in AGN (Perna et al. 2017; Brinchmann et al. 2023). This is 
further illustrated in Fig. 3, where we plot nearby AGN (red contours) 
and star-forming galaxies (grey contours) from SDSS (Abazajian 
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Figure 3. Narrow line ratio diagnostic diagram of [O 11]4.4363/H y versus 
[O 11]45007/[O m]à4363. Different regions within ZS7 are indicated by 
golden symbols: the BLR location (star), the [O 111]A5007 centroid (diamond), 
and the surrounding region (circle; see Section 7 for details). Local star- 
forming galaxies (grey shading and points) and AGN (red shading and points), 
dust-corrected following Cardelli, Clayton & Mathis (1989), separate into 
two parallel sequences in this diagram, with AGN occupying a region of 
higher electron temperatures and having elevated ratios of [O m]à4363/H y. 
The red diamond shows a galaxy from the JWST early release observations 
(Pontoppidan et al. 2022) at similar redshift, which is identified as an AGN 
through a [Ne 1V]A4.2422, 2424 detection (Brinchmann 2023; B23), and the 
red star indicates the narrow line ratios of a BLR AGN at z — 5.55 (Übler 
et al. 2023; 023). The BLR centroid of ZS7 is located at the extreme end of 
the local AGN distribution, and also the [O 11]4.5007 centroid falls on top of 
the local AGN cloud, while the host emission is located in a region where 
local star-forming galaxies and AGN overlap. The narrow line ratios extracted 
from the [O 11]A5007 centroid position could indicate recent AGN activity 
also in this part of ZS7. 


et al. 2009) in the [O 111]44363/H y versus [O 111]45007/[O 111114363 
line ratio diagram. In this diagram, local star-forming galaxies and 
AGN separate into two parallel sequences, with AGN occupying a 
region of higher electron temperatures and having elevated ratios 
of [O 11]44363/H y (Fig. 3; also Mazzolari et al., in preparation). 
Indeed, the narrow emission line ratios at the location of the broad 
H f in ZS7 lie at the extreme end of the nearby AGN distribution 
(golden star in Fig. 3). This is yet another confirmation that the broad 
H f is tracing a (Type 1) AGN. 

Intriguingly, as illustrated by the bottom-right panel of Fig. 2, the 
auroral [O 111]A4363 line is spatially extended across the locations 
of both the BLR and [O n]à5007 centroid. We can measure the 
[O m]à4363/H y and [O m]à5007/[0 11]A4363 line ratios at the 
[O 11]45007 centroid position, and we find that also the emission line 
properties at the location of the [O 11]A5007 centroid are consistent 
with local AGN, as shown by the golden diamond in Fig. 3. This 
finding strongly suggest that also the galaxy associated with the 
[O m]à 5007 centroid is hosting a (Type 2) AGN. 

We finally note that, if we extract these emission line ratios from an 
integrated spectrum excluding the region where the BLR is detected 
with S/N > 5, and excluding the 3 x 3 spaxels around the [O 111]45007 
centroid, then we find values falling into the region which is locally 
occupied by both AGN and star-forming galaxies (golden circle in 
Fig. 3). 
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8 AN IMMINENT GALAXY AND MASSIVE 
BLACK HOLE MERGER 


As discussed in the previous sections, the NIRSpec data and the 
NIRCam data provide evidence that ZS7 is composed of two 
components in close vicinity. The projected separation of 620 pc 
and velocity difference of the narrow line emission between the BLR 
location and the centroid of the [O 11]A5007 emitter (720 kms“), 
likely hosting another AGN, are small, and likely imply an imminent 
merger. 

In Section 5, we measured the mass of the black hole associ- 
ated with the Type 1 AGN, traced by the broad H f, by using 
virial relations. Estimating the mass of the putative second black 
hole, associated with the [OrJ]A5007 emitter, is more difficult. 
Using our stellar mass estimate of the [O m]à5007 emitter (M, 
^ 2.9 x 10? Me, ie. excluding the region where the BLR is 
significantly detected, see Section 3.3) and assuming a M,—M, 
relation calibrated on nearby AGN (Reines & Volonteri 2015), results 
in a secondary black hole mass associated with the [O m]à 5007 
emitter of log(M, /M) ~ 5.8 + 0.6. If we assume a lower limit for 
the stellar mass of M, ^ 1.7 x 10? Mo (see discussion in Section 3.3), 
we find log(M./Mo) ~ 5.6 + 0.6. Comparing this to the black hole 
mass estimate for the BLR AGN (log (M,/Mo) = 7.7 + 0.4), we 
find a ratio that is very similar to the dual AGN candidates identified 
by Maiolino et al. (2023). If instead we assume a M,—M, relation 
calibrated on local S/SO and elliptical galaxies, we get a black holes 
mass higher by about one order of magnitude. Alternatively, if we 
assume a black hole-to-stellar mass ratio between 0.01 and 0.1, 
as found for many high-z BLR-AGN discovered with JWST, the 
corresponding mass of the putative second black hole would be 
comparable to or even larger than the mass of the BLR-identified 
AGN in ZS7. We note that simulations find that in ‘offset? AGN 
systems (i.e. interacting galaxies in which only one of the two black 
holes is accreting) the active black hole is typically more massive 
(Steinborn et al. 2016; Chen et al. 2023), although these results are 
only available at z < 3. 

Continuing with the interpretation of an imminent massive black 
hole merger, at a projected separation of 620pc the black hole 
pair is likely in the dynamical friction phase when the interaction 
with stars and gas assists the black holes in reducing their angular 
momenta to sink to the common centre of mass (Chandrasekhar 1943; 
Ostriker 1999), before three-body interactions with individual stars 
(Begelman, Blandford & Rees 1980; Quinlan 1996) and eventually 
the emission of gravitational waves lead to coalescence of the black 
hole pair (Peters & Mathews 1963; Peters 1964). Using our black 
hole mass estimate, the stellar mass estimate of the [O m]à5007 
emitter summed within roughly 3 kpc, and assuming an absolute 
separation of the BLR and the [O 11145007 centroid of 1 kpc, we find 
an approximative inspiral time of 100-200 Myr (Binney & Tremaine 
2008). 


9 IMPLICATIONS 


It is not possible to provide any reliable statistics on the number and 
fraction of merging black holes at high redshift, based on this single 
case. We note, however, that this galaxy was observed as part of a 
sample of 11 galaxies at z > 6 that were observed with the NIRSpec- 
IFU mode within the GA-NIFS survey. Of these, the ZS7 system and 
a second galaxy (Cresci et al., in preparation) are found to clearly 
host an AGN based on the detection of broad Balmer lines. The other 
AGN at z > 6in GA-NIFS is also found in an interacting system, on a 
few kpc separation. Also the two z ~ 6.8 QSOs analysed so far within 
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GA-NIFS undergo galaxy-galaxy mergers (Marshall et al. 2023; see 
also e.g. Yue et al. 2021). Although these are low-number statistics, 
these results suggest that AGN (and more broadly black hole activity) 
in the early universe are associated with galaxy interactions. In the 
case of ZS7, the separation of the individual systems is so small 
that the galaxies and their supermassive black holes will most likely 
merge within the next few hundred million years. At slightly lower 
redshifts (4 — z < 6), a recent study by Maiolino et al. (2023) 
identified a significant number (3/11) of candidate dual AGN based 
on the complex profile of broad Balmer lines in the NIRSpec multi- 
object JWST Deep Extragalactic Survey (JADES; Eisenstein et al. 
2023, without 2d spatial information). At z ~ 3, Perna et al. (2023b) 
identify 20-30 per cent dual AGN on separation scales of 3-30 kpc 
within GA-NIFS. Summarizing, the findings of our and other studies 
using different approaches to assess black hole activity in the early 
universe are consistent with a potentially large number of black holes 
in the process of merging. 

There are no simulations available that provide the expected 
number of binary black holes, with separations less than | kpc, in 
this mass range and at z > 7. However, the mentioned finding of a 
large fraction of dual AGN at lower redshift (z~3) on scales of a 
few/several kpc already highlighted tension with the expectations of 
simulations, which predict a lower fraction of dual AGN at z « 5 (De 
Rosa et al. 2019; Rosas-Guevara et al. 2019; Volonteri et al. 2022; 
Chen et al. 2023) than observed (Koss et al. 2012; Spingola et al. 
2019; Perna et al. 2023b). However, while some simulations find 
an increasing fraction of massive black hole pairs with decreasing 
redshift due to long coalescence time-scales of massive black hole 
mergers in these models (Chen et al. 2023), most simulations predict 
the fraction of massive black hole pairs to increase with increasing 
redshift due to higher galaxy-galaxy merger rates (Volonteri et al. 
2016, 2022; De Rosa et al. 2019; Rosas-Guevara et al. 2019; see also 
discussion by Sesana, Volonteri & Haardt 2009). 

We also mention another possible scenario explaining our observa- 
tions of ZS7, namely that the displaced accreting black hole actually 
results from the recoil of a recent supermassive black hole merger. 
Anisotropic emission of gravitational waves during a black hole 
merger event may give the resulting black hole a velocity kick (Peres 
1962; Bekenstein 1973). Some simulations predict these velocities 
could be as high as several hundred to a few thousand km/s for pre- 
merger black holes with comparable masses, which could even lead 
to the ejection of the merged black hole from the host galaxy (e.g. 
Madau & Quataert 2004; Merritt et al. 2004; Blecha et al. 2016). 
Recoiling black holes may be observable through both a spatial and 
a velocity offset to the host galaxy (see e.g. Castignani et al. 2022), 
and this could be the case in ZS7: a spatial offset of 620 pc between 
the host galaxy emission and the BLR is evident in this case. In 
addition, we register a small velocity offset between the narrow line 
emission of ZS7 and the BLR emission of ~40 km s^! . While recoil 
velocities depend on the masses and relative spins of the initial black 
holes, the fact that we do not observe a large velocity offset of 
the BLR makes the interpretation of a recoiling supermassive black 
hole less appealing, although we caution that we cannot measure 
the absolute velocity difference. Furthermore, a recoiled black hole 
would probably not have a NLR around it and, if it did, it would 
likely be the interstellar medium of the host galaxy, therefore it 
should share the same velocity field as the host. In the case of ZS7 
we clearly see that the narrow component of the nebular emission 
around the accreting black hole, its NLR, has a different velocity 
(Fig. 2). This suggests a more likely scenario in which the accreting 
black hole has yet to merge and is still embedded in the interstellar 
medium of its own galaxy. Finally, the indication that the [O 111]4.5007 
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centroid is also likely to host another (Type 2) AGN, implies that the 
merging has yet to happen, unless the narrow lines signatures are 
actually the fossil record of past accretion at that location. However, 
even the recoil scenario would be pointing at a black hole merging 
event; the difference between the two scenarios would be that in one 
case the merging has yet to happen, while in the other case it has just 
happened. 

Finally, we note that another possible explanation for the offset 
BLR could be that it is tracing one of the so-called ‘wander- 
ing IMBHs'. Recent simulations have shown that (at least on 
scales of several kpc) the dynamical friction during galaxy inter- 
actions/mergers may not be efficient at sinking black holes at high 
redshift (Di Matteo et al. 2023b). As a consequence, this may result 
into a potentially large population of IMBHs in the haloes of galaxies. 
However, such a population of wandering IMBHs is expected to be 
on scales much larger (tens of kpc) than the separation observed in 
ZS7 (where dynamical friction is much more efficient). Additionally, 
a ‘wandering’ black hole, even if accreting, should be essentially 
‘naked’, that is, should not have its own NLR. This is in contrast 
with what is observed for the offset BLR in ZS7, which is clearly 
exciting at least [O 11]A4363 tracing low-density gas outside the 
BLR. 


10 CONCLUSIONS 


The results presented in this work highlight the power of NIRSpec- 
IFU observations in the study of high-z galaxies, mergers, and 
the detection of moderately massive black holes through imaging- 
spectroscopy. Our observations provide clear and robust evidence 
for a massive black hole involved in a merger with another galaxy, 
likely hosting another accreting black hole, at z = 7.15, only 740 Myr 
after the big bang. Overall, our results seem to support a scenario 
of an imminent massive black hole merger in the early universe, 
highlighting this as an additional important channel for the early 
growth of black holes. Together with other recent findings in the 
literature, this suggests that massive black hole merging in the distant 
universe is common. Our observations may be used as guidance for 
the modelling of gravitational-wave events originating from massive 
black hole mergers that will be detectable with future observatories 
like LISA (e.g. Haehnelt 1994; Jaffe & Backer 2003; Sesana et al. 
2005; Valiante et al. 2021; Amaro-Seoane et al. 2023). 
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APPENDIX A: FULL SPECTRAL FITS 


In Fig. Al we show full spectral fits to spectra extracted from 
individual spaxels at the BLR location, and at the centroid location of 
the [O m]à5007 emitter. At the BLR location, the broad H 6 emission 
is evident. Some flux of the BLR point-source is spread through the 
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Figure A1. Full spectral fits to spectra extracted from individual spaxels at the BLR location (top), and at the centroid location of the [O 11]A5007 emitter 
(bottom). Residuals show (data — model)/uncertainties. We show our full fit in blue, with narrow emission lines in green and the BLR component in purple. At 
the BLR location, the broad H £ emission is evident, while the [O 111]AA4959, 5007 emission is narrow at all locations. We note the strong auroral [O 111]A4363 


emission line in both spectra. 


PSF (PSFrewum ~ 0.13 — 0.14" at 4 um) and therefore is still visible 
at the [O 111]A5007 centroid position. 


APPENDIX B: DUST EXTINCTION 


In Fig. B1 we show a map of the PSF-matched H y to H 8 narrow 
emission line flux, derived using the empirical fitting functions for 
the NIRSpec-IFU PSF provided by D'Eugenio et al. (2023). The 
centroid locations of the [O 11].5007 and BLR are indicated by a 
cross and star, respectively. The black line in the colourbar shows 
the theoretical line ratio of H y/H 8 — 0.468, whereas lower values 
indicate dust extinction. The data show that there is no enhanced 
dust extinction at the location of the BLR. This rules out a scenario 
where the BLR centroid location is the true, dust-obscured centre of 
the ZS7 system. 


€ 2024 The Author(s). 


D ow) uds 


1 kpc 


EE i, 


Figure Bl. Map of the PSF-matched H y to H f narrow emission line flux 
where both lines have S/N > 5, with the centroid location of the [O 11]A5007 
and the location of the BLR indicated by a cross and star, respectively. 
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